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,—1 . Abstract: We calculate the thermal Euclidean correlators and the spectral functions of the energy- 

^ momentum tensor for pure gauge theories, including at non-zero spatial momentum, at leading order 

I in perturbation theory. Our goal is to improve the extraction of transport properties from Euclidean 

^ correlators that are computable in lattice QCD. Based on our results and the predictions of hydrody- 
namics for the structure of the spectral functions at low frequencies, we show that the shear and bulk 

Q viscosities can advantageously be extracted from the Euclidean correlators of the conserved charges, 

■ energy and momentum, at small but non-vanishing spatial momentum. The spectral functions in these 

• • . channels are free of the ultraviolet u;^ term which represents a large background to the thermal physics 

• , encoded in the correlators of the fluxes. 
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1. Introduction 

The particles produced in RHIC heavy-ion coUisions exhibit sizeable elliptic and radial flow [1, 2, 3, 4]. 
This is a signature of collective behavior, which the equations of ideal hydrodynamics have been able 
to describe quite successfully [5, 6, 7]. The leading corrections in a gradient expansion of flow velocities 
is parametrized by shear and bulk viscosity, rj and Detailed viscous relativistic hydrodynamics 
calculations [8, 9, 10, 11, 12, 13] find that the experimental data excludes a shear viscosity to entropy 
density ratio larger than about 0.3, in agreement with early rough estimates [14]. This would make the 
substance created an exceptionally good fluid [15], and it is therefore of primary interest to improve 
our knowledge of the transport coefficients of the quark-gluon plasma. 

Theoretically, the transport coefficients can be extracted from the imaginary part of a retarded 
two-point correlation function; this imaginary part is called the spectral function and is denoted by 
p{(jj). This is the content of the Kubo formulas [16]. By analytic continuation, the spectral function 
is also related to Euclidean correlators Ce [17, 18], which are functions of Euclidean time xq, spatial 
momentum q and temperature T, by 

CE(xo,q,T)= rdup{u,ci,T) ^"^^^(^0 -2t) _ (11) 
Jo smh 

The spectral functions are odd in uj and have everywhere the same sign as w. Since viscosities 
parametrize the dissipation of momentum, the relevant operators are elements of the energy-momentum 
tensor T^iy. Because the temperatures reached at RHIC are not much larger than twice the critical 
temperature Tc where a gas of hadrons rapidly crosses over to a system with many more degrees of 
freedom (see [19] and Ref. therein), perturbative methods [20, 21, 22, 23, 24] are not directly applicable. 
They can however tell us about the asymptotic high T behavior of the viscosities and the associated 
spectral functions. They also allow us to familiarize ourselves with the intricate kinematics that arise 
at non-zero temperature and spatial momentum, one of the objectives of this work. Finite-momentum 
meson spectral functions were investigated analytically in [25]. 

Lattice calculations of viscosities [26, 27, 28, 17] have so far mostly focused on the Euclidean 
correlators of T12 for the shear viscosity, and Ta (or T^^) for the bulk viscosity. Indeed the Kubo 
formulas read in these cases 

ry(r) = vr lim , (1.2) 

t<;-+0 U! 

C(T) = f lim /"'"'"-"-^' , (1.3) 

where Pfj,u,pa corresponds to {Tfj^uTpa)- As we shall review in section 3, hydrodynamics, as an effective 
theory describing low-frequency phenomena around equilibrium, predicts the functional form of the 
spectral functions in these channels, including at small but non-vanishing q. 

Leaving aside the problem of determining the spectral function given the Euclidean correlator, a 
particular difficulty in this approach is that the spectral functions grow as u'^ (times a power series 
in as) at large frequencies. This buries the contribution from the small lo region under a much larger 
contribution (by a factor of at least five) from ultraviolet modes. The latter contribution is almost 
temperature independent and therefore does not advance in any way our understanding of thermal 
physics. Given that the Euclidean correlators are determined by Monte-Carlo methods and carry 
statistical errors, this is a significant drawback. The difficulty is far more severe than in studies of the 



charmonium spectral functions [29, 30, 31], or calculations of the electromagnetic conductivity [32], 
because the spectral function only grows as uP' for the vector current V'7iV'- 

Methods to subtract the ultraviolet contribution to the Euclidean correlator, and hence to enhance 
the sensitivity of the lattice observables to the low-frequency region described by hydrodynamics, have 
been proposed and implemented [18]. One of these methods (subtracting the T = spectral function) 
has the virtue of removing the contribution completely, by contrast with a perturbative order-by- 
order subtraction. The drawback is that positivity of the integrand in Eq. 1.1 is given up, and a large 
part of the signal is lost in the difference, an unfavorable situation from the numerical point of view. 

Here, based on the exact Ward identities that follow from the conservation of the energy-momentum 
tensor, we show that the spectral function of the energy density operator with non-vanishing spatial 
momentum goes to a constant as ^ oo. Similarly, the two-point function of the momentum density 
operator grows only as llP' . This is confirmed by our perturbative calculation. These Euclidean cor- 
relators are therefore far more sensitive to the thermal effects than the correlators of the fluxes. Yet 
provided q is sufficiently small, their low-frequency region is still described by hydrodynamics, and 
therefore the shear and bulk viscosity can be extracted from them. 

We collect leading perturbative results on the Euclidean correlators and spectral functions in 
section 2, and discuss in section 3 the interplay of the perturbative predictions with the hydrodynamics 
predictions to propose a new way to extract shear and bulk viscosity from Euclidean correlators. 



2. Perturbative calculation 

The Euclidean energy-momentum tensor for SU(A^c) gauge theories reads 



T^.(.x) = + Jv^(a^) (2-1) 

a I \ _ 1 r pa pa _ pa po 9) 

^(^) = ^ (2-3) 

= -M' + ---, &o = Sk. (2.4) 



In the U(l) case, the summation over the adjoint index a is trivial and of course jii^g) = 0. In contrast 
with Minkovsky space, Toj = % is an antihermitian operator. In particular, (roi(2;)Toi(0)) < for 
X 7^ 0, and Pj = i J (fyLTQj{x) is the usual momentum operator, for instance -P, ]q) = 9jlq) for a 
one-particle state. 

In view of the form of the energy-momentum tensor, we consider now the connected correlators of 
the field strength tensor and those of the field strength tensor squared, at leading order in perturbation 
theory. The expectation values obtained in that approximation are denoted by (. . .)o and throughout 
this paper we keep only connected diagrams. The effect on spectral functions of using Hard-Thermal 
Loop resummed perturbation theory [33] in the region a; < T is beyond the scope of this work. 
Nevertheless our results are perfectly physical, since they are exact in the U(l) case. 

In Feynman gauge, we have 

{F^„{x)Faf3{y))o = dA (t)^ua0{x - y) , (2.5) 

(I>nuaf3{x) = 5vpfna{x) + S/^afupix) - 6^pf^a{x) - 6^affip{x) , (2.6) 

U^{x) ^Tj^J (2^)3 p2 + p2 P-P^- (2-7) 



One then easily obtains 



{F'llAx)Fpaix) F!^p{y)F!^g{y))o = dA 4>^uap{x - y)4>p<j^5ix - y) + (l)^u-rd{x - y)4>paa(}{x - y) ■ (2.8) 

The color factor is (Ia = N'^ — 1 for SU (iVc) and 1 for U(l). To study mixed correlators (which 
are functions of (xo,p)), we introduce the spatial Fourier transform of (f)fj,i>ai3ix), (j^nvapiPiXQ) = 
I d^^cf>^,ap{x) e^P ". Then 



(2.9) 



= dA 



P 



(pfiiyapiP, Xo)<j)pa^s{-{p + q), Xq) + (pf^y^siP, Xo)(t>paap{-{P + Q), a^o) 



The integral is worked out in appendix A. The results are of course gauge-invariant. In the next 
two subsections, wc present the results in a number of channels of physical interest. Throughout this 
section, we align the spatial momentum with the z-axis, q = qe^, q > 0. Since the spectral function 
is odd in uj, we also choose a; > without loss of generality. The most general results are given at 
the end of the section, however we find it useful to give the explicit form in the simpler cases of zero 
temperature or zero momentum. 



2.1 Zero temperature 

At T = and zero spatial momentum, by dimensional analysis, the correlation functions at tree-level 
all fall off as l/xg. With finite-momentum and for xq > 0, 



I 



d^xe^^-^(ri3(x)ri3(0))o 

Pl3,13(<^,g) 

d^^e'^-^{Tuix)Tumo 

Pi2,i2('^,g) 
J d^^e'^-^{e3s{x)e3s{0))o 

P33,33('^,g) 

j d=^xe^'i-^(0(x)0(O))o 
J d^^e''^-^{eoo{x)9oomo 



dAe~ 



qxo 



{qxQ + 2) {q Xq + SqxQ + 6) , 



5(47r)2x^ 
dAO{uj-q) 2/ 2 2x 

10(4.)^ 
4dAe-«^o , 2 2n 

.f. .2 5 ^ + qXo + q Xq) 

dAOjuJ - q) ^ 2 „2^2 
10(47r)2 

2^4 6-9^0 



/ 2 2\2 

-(w -q) , 



5 (24 + 2^XQ + 12q^xl + 4g=^xg + q^x^), 



15(47r)2xo 
2dAe{u - q) 4 
15(47r)2 ' 

UasN.V 2dAe 



Git 



-qxo 



(47r)2x3 



llasNc\ dAO{ijJ - q) , 2 2 



67r J 4(47r)2 
2(i^g^e-«=^o 
15(47r)2xo 
2dAe{u - q) 4 
15(47r)2 ^ ■ 



(3 + 2,qxo + q^xl) , 



(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 



2.2 Finite temperature 

The q = correlation functions of interest are 



J d^x (eoo(0)eoo(a;))o = 
J (fx (^(O)0oo(a^))o = 0, 
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\ 67r 

32dAT^ 



J d^x {e{o)e{x))o 

J d^X (^12 (0)^12)0 

/<,.(M0).n)„ = i^(/«-5) 



where r = 1 — 2Txc] and 



/(r) = / 

^0 



°° , 4 COsh^ TS TT^ 3 



n 



+ 



n 



sinh^ s 60 4 4-; (n - r)^ (n + r 

n>l 



(2.20) 
(2.21) 

(2.22) 
(2.23) 
(2.24) 

(2.25) 



The two point function of J" (i^xroo(x) is time-independent, as expected for a conserved charge. The 
spectral functions in the tensor and scalar channels read (see for instance [34, 26]) 

fl A I, A / 9-7r\ ^ 

Pi2,i2('^,o,r) 



2.2.1 Finite momentum 

For a polynomial P, we define 

I{[P],u,q,T)=e{u-q) C dz 

Jo 

Then the spectral functions read 
Pi3,i3('^,g,r 

Pi2,i2(w, r 

P33,33('^,g,r 
P0,e('^,g,r 

poo,oo(w,?,r 
-poi,oi('^,g,r 

-P03,03(<:^,g,r 



10(47r)2 tanh ^ V 15 
dA ( \\a.sNf\ ^ u 



4(47r)2 V 67r y tanh^' 



(2.26) 
(2.27) 



P{z) sinh 



2T 



cosh ^ — cosh ^ 



+ 



°^ P(z) sinh 
dz 



2T 



cosh ^ — cosh ^ 



8 (47r)2 

dA 



32 (47r) 



U;2(a;2-g2) X([i_^4j^^^^^2^^^ 
J {iv^-q^f I{[l + 6z^ + z%iV,q,T), 



dA 



4(47r) 



X([(l-zY],a;,9,r), 



(u;2-g^)^ Ii[l],uj,q,T), 



4(47r) 



67r ; 4(47r)2 

2 J([(l-z2)2],u;,(^,T), 



X([(l-z2)2] 

g2(a;2-g2) J([(l - /)], a;, T) , 



8(47r: 



i([(i-.Y],^,9,r). 



4(47r) 



(2.28) 

(2.29) 
(2.30) 
(2.31) 

(2.32) 
(2.33) 
(2.34) 
(2.35) 




One finds that 

X l,a;,g,r =— %-a; + _ log --^ 2.36 

q q smn |a; — (/|/4J 

and hence the trace anomaly correlator can be expressed in terms of elementary functions. The integrals 
with and in the numerator can be expressed in terms of polylogarithms, explicit formulas are 
given in appendix C. A few spectral functions are displayed on Figure 1. 

3. Physics discussion 

At small momentum and frequency, the expression for the spectral functions of the momentum densities 
are predicted by hydrodynamics (see [35] for an explicit derivation), 



(3.1) 
(3.2) 



where s is the entropy density, Vg is the velocity of sound and q = qe^,. 

Based on the fact the matrix elements of d^T^^^ vanish between any two on-shell states, the 
Euclidean correlators of the charges and those of the fluxes are related in a simple way. We emphasize 
that these relations are exact, since they derive from a Ward identity. In terms of the spectral functions 
they read (q = qez) 

<^^Poo,oo('^,q) = 9^33,33 (<^,q) (3.3) 

-t^^poi,oi('^,q) = g^pi3,i3('^,q) (3.4) 

-t^^P03,03('^,q) = g^P33,33('^,q). (3.5) 



These relations are in particular satisfied by our trcclcvcl expressions, Eq. (2.29-2.35). Note that the 
minus signs on the right-hand side of Eq. (3.4-3.5) are absent in Minkovsky space (they come from the 
definition of Tofe itself, see the remark at the beginning of section 2). 

Equations (3.4-3.5) and (3.1-3.2) can be combined to obtain (1.2) and (1.3), which have so far 
been the basis of the calculation of shear and bulk viscosity using lattice Monte-Carlo techniques [26, 
27, 28, 17]. In that strategy, the momentum q is set to zero at the outset, and u is sent to zero at the 
end. 

However, in view of Eq. (3.1-3.2), the shear and bulk viscosity can be extracted from the low- 
freqency behavior of the spectral functions for the four charge densities Tq^, as long as q 7^ 0. The 
advantage of using these correlators is that the ultraviolet contributions are highly suppressed compared 
to the correlators of the spatial components. We may illustrate this point by two numerical examples. 

In Eq. 2.26, relevant to shear viscosity, the w^/tanh ^ term contributes for 86% to the Euclidean 
correlator at t = 1/2T. By contrast, in the /9oi,oi channel for q = '7rT/2, the contributions to C-e{x() = 
1/2T) / {(IaT^) coming from uj > q and to < q are respectively ^ 0.04 and 0.6, assuming the treelevel 
form (Eq. 2.34) in the first region and Eq. 3.1 with s = |ssb = jgdAT^'^T^ and 77/s = l/47r in the 
second. 

In the energy density channel, the increase in sensitivity to the low-frequency region of the spectral 
function is even more dramatic: for q = vrT/2, the contribution to Ce{xo = \/2T) / {(IaT^) from uj > q 
is merely ~ 0.01 based on Eq. 2.33, while the sub-threshold contribution is about 1.9 for s = l^sBj 
v1 = i, rj/s = l/Air and ( = 0. These values are inspired by the strongly coupled Af = A SU(iVc) 
SYM gauge theory, which can be studied by analytic AdS/CFT methods (see for instance [36] and Ref. 
therein). 

One of the key issues in practice is to achieve sufficiently small q for the hydrodynamics prediction 
to be valid below threshold. One will want to reach q < ttT/2 and check this explicitly. Since q = 2tt/L 
is the smallest momentum available in a periodic box, this requires simulating in rather large spatial 
volumes. Approaching a second order phase transition, it may become impractical to reach sufficiently 
small momenta. Anisotropic lattices [37, 18] can help to achieve large physical volumes while keeping 
discretization errors under control. Prom the algorithmic point of view, having a non-zero spatial 
momentum is a particularly favorable situation for the multi-level algorithm [38, 39, 40], since it allows 
the ultraviolet fluctuations to be tamed very efficiently. Indeed non-perturbative, non-zero momentum 
correlators of the momentum fluxes were presented with 1% precision in [18]. 

In the case of the vector current tp^^ip, it may also be profitable to exploit the correlator of the 
charge density with non-zero momentum. Non-zero momentum correlators of the current have been 
computed with good precision [41], exploiting twisted boundary conditions to scan low momenta more 
easily. 
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A. Calculation of Euclidean correlators 

In this appendix we derive a form of the Euclidean correlators which is well suited for accurate numerical 
evaluation. We start with (ri2(0) ri2(x)), and then give results for other channels treated in the same 



way. The orientation of the momentum is specified by q = qei, i = 1,2,3. The color factor dA is 
omitted. 

A.l The (ri2(0)ri2(x)) correlator 

(ri2(0)ri2(x)) = ((i=^10i^20)(0) {FioF2o){x)) + ((Fi3F23)(0) (Fi3F23)(x)) 

+2 ((Fl0i^20)(0) {Fi3F23){x)). (A.l) 

Exphcitly, the tensor 4>fxval3 reads 

In the r — > Hmit of course one can make the substitution ~^ I Defining 

— 1 .r-^ /" d^p g^ipo^o gigolo 



we have 



y" d=^xe^'i-"(FioF2o)(0) (Fioi^2o)(x)) = ^ + B , (A.6) 

y d3xe*^"(Fi3i^23)(0) (Fi3F23)(x)) =B + C, (A.7) 

y" d=^xe^'i-"(FioF2o)(0) (Fi3i^23)(x)) = P . (A.8) 
Using the Poisson summation formula, we obtain 



A = 


m,n6Z 


B = 


m,n6Z 


C = 


is/ 

m,n6Z 


P = 


m,n£Z 



^ Vi^l - ; U/^2 ^ ^^2; - Vt^ ^ ) -|p||xo+mLo| -|p+q||xo+nLo| fA Q) 

(27r)3 IpI |p + q| ' V • ; 

^^P yiP2(yi+gl)(p2+g2) ^-|p||a;o+mLo| g-|p+q|No+nLo| 10) 

(27r)3 |p||p + q| ' 



y ^P3(P3 + g3) e-|Pll«^ole-|P+'»ll-°+"^°l, (A.12) 
where 

s„(x) = sign(x + n) . (A. 13) 



Integrating in spherical coordinates leads to 

^5 /-oo 



where 



m,neZ ^ 

X + 1 (J4 + 2 + Jc) [q\xQ + nLol Vfc^ + 1, (A.14) 

^Afi,c(x,AA;) = r ^^e-^^PA,B,c{ylA,k) (A.15) 

/•^ , 

iT,(x,AA;)= J ^dye-'^^Pviy/Ak). (A.16) 



Via the change of variables u = s/l + y, these integrals can be obtained in terms of the function 

g{x, A) = i(e-^^T^ - e"^^^). (A.17) 

We have 

JA,B,c{^,A,k) = 2Pa,b,c{a-\£, - g{x,A) (A.18) 

Iv{x,A,k) = -2^Pv[A-\^-l),k)gix,A) (A.19) 

Note that ^/T^ = J^$iL for A = tJtt- 

One then obtains the correlator in the form 



,5 

n,meZ 



j d^xe'^-^{Ti2{0)Ti2{x)) = 5 ^ g ([n]; s„,(f^) s„(f^), g|xo + mLo|, q\xo + nLo\) , (A.20) 
where 



g([n];.,x,y) = ^-'^('^yj^ e-J^n(.,y,x) ^ g([n];.,y,x) (A.21) 

\x y ) 

and the appropriate polynomials 11 will be specified below for a few cases of interest. The function 
g will always be finite at x = y, and we give the first two terms in the expansion around this point, 
which is useful when implementing expression (A.20) numerically. Since the series in n and m are 
exponentially convergent, only a few terms are needed to reach, say, a precision of 10~^. 

Consider the case q = (g, 0, 0), q> 0. Relabelling the p-coordinates pi — ps, Pa — P2 , P2 — Pi in 
Eq. A. 12 before going to spherical coordinates, we obtain, 

PA{y, k) = il- 2/2) (^y^ + ^y+^ + l^^ (A.22) 

PB(y,A;) = y(l-y2)(y + ^), (A.23) 

Pc{y,k) = l-y\ (A.24) 
Pv{y,k) = l-y\ (A.25) 



and from there 

2nf^i2(s>^>2/) = 4608 + x(6 + a;)(768 + a;(240 + x(72 + a;(10 + x)))) (A.26) 
-4sx(48 + x(48 + x{12 + x)))y 
-(96 + x(96 + x{120 + x(32 + 3x))))y'^ 
+8sx{6 + x)y^ + (6 + x)(-2 + 3ar)y^ - 4sy^ - / , 

g([n^^^i2]; s,x + le,x- ie) = ^££^|x3(2 + s) + 3(7 + 5s) + 3x(7 + 5s) + x\9 + 6s) 

2 

^ ((39 + 21s + x^{A + s) + 3x(13 + 7s) + x^il7 + 8s))) + 0(e^)} . 

In Eq. A. 27, we have defined g by continuity. 
For the case q = (0, 0, g), g > 0, we obtain 

„.,x l^So/l 3\2 2 1 1 

Prs{y,k) = l{l-y')\ 
Pc{y,k) = PA{y), 
Pv{y,k) = 2y{y + ^), 

and further 

Ull^^^{s,x,y) = -576-x(4 + a;)(144 + x(48 + x(24 + x(6 + x)))) 
+2sx(96 + x(96 + x(48 + x(ll + x))))y 
+(-48 + x(-48 + x(-48 + (-4 + x)x)))y^ 
-4sx^(5 + x)y^ + x(14 + x)/ + 2s(-l + x)y^ - y 

g{[^%,i2\^s, X, x) = 2^ {21 + x(21 + x(9 + x(2 + x))) 
+s(15 - x(-15 + (-1 + x)x(3 + x))) 

— ( 69 - 2x^(-l + s) - x^(-l + s) + 63s 
56 V 

+x2(25 + 19s) + x(69 + 63s)) + 0(e'')| . 

Note that the polynomials PA,B,C,v{y, k) all have the symmetry P4,e,c,x>(-y, -k) = PA,B,C,v{y, k). 

A. 2 (TiiTii) with momentum along x direction 

For / d3xe*'i ''(rii(0)rii(x)), the polynomial is 

-Ul\n{s, X, y) = 4608 + x(4608 + x(2112 + x(576 + x(108 + x(14 + x))))) 
+sx(2 + x)(120 + x(60 + x(12 + x)))y 

-(192 + x(192 + x(24 + x(4 + x))))y^ - 2sx(18 + x(8 + x))y^ 

-(-12 + x(10 + x))/ + s(2 + x)/ + / , (A.34) 



(A.27) 



(A.28) 

(A.29) 
(A.30) 
(A.31) 

(A.32) 



(A.33) 



g([nfl,n];s,x + ie,x- |e) = (A.35) 



|66(1 + x) + 2a;2(15 + x(4 + x)) + s(30 + x(30 + x(15 + x(5 + x)))) 

2 

+|g (^138 + 138x + x^(2 + s) + x2(8 + 5s) + x^(54 + 5s)) + 0(e^)} . 

A. 3 Scalar correlator with non-zero momentum 

For iE;.<.,p<a/^'xe^'''"(^M.(0)^pa(^)), t^e polynomial is 

-npg{s, X, y) = x^(12 + x(6 + x)) - 2sx^(24 + x(9 + x))y (A.36) 
-x^(-72 + (-12 + x)x)y2 + 4sx(-12 + .x(3 + x))y^ 
-(-12 + x(18 + x))y^ - 2s(-3 + x)/ + 

g([n^^,]; s, X + ^e, X - \e) = { - x\-l + s) + 45(1 + s) + 45x(l + s) + 15x^(1 + s) 

( - x^(-l + s) + 105(1 + s) + 105x(l + s) + 35x^(1 + s)) + 0(e^)} . (A.37) 

B. Calculation of spectral functions 

In this appendix, we show how to obtain Eq. (2.29-2.35) in the case of the (T12T12) correlator. Other 
cases can be treated in exactly the same way. In this section q is understood to be in units of 2T. 
Using the propagator in the mixed representation, 

1 ^iPQXQ ^ 1 cosh|p|(iLo-xo) 

Lo^p^o + P^ 2|p| sinhi|p|Lo ' ^ ' ' 

we obtain (r = 1 — 2Txq) 

A,B,C = HpUp C dx / Pa,b,c{x,^) (B.2) 

tt^LIJq J_i p^ + q^ + 2pqx ' ' ^ 1' 

cosh(p + \J p^ -\- cp- -\- 2pqx)T + cosh(p — \/p^ + q^ + 2pqx)T 
cosh(p + \/p^ + q^ + 2pqx) — cosh(p — \/p^ + q^ + 2pqx) 



For the term proportional to cosh(p + ^/p"^ + q^ + 2pqx)T, we do the change of variables p = 2{k+qx) ' 
which brings this expression to the desired form cosh(A;r) on the integration interval k\'^ . 

For the cosh(p — ■\/ p^ + q^ + 2pgx)r term, we choose p = 2{qx-k) ' 'w^^ich achieves the same on the 
integration interval k\qx- We then notice that if J^{—x, —k, q) = J^{x, k, q), 

/I f-q i^q 

dx / dkJ^{x,k,q) = / J^{x,—k,q). 
-1 Jqx Jo J-i 

One therefore finds that the first term determines the spectral function above the threshold q, while 
the second determines it below q. We write 

p{lo, q, T) =p<{u, q, T) + p>{io, q, T), (B.3) 



where the first term vanishes for a; > g and the second for a; < 



After some algebra, 



PA,B,C,> i^kT, 2qT, T) = ^ ^^^^ ^ ^ (B.4) 



{k + qxY coshA;-coshg|±g- 

and after the further change of variables x = we obtain 

^ „ ^ ^, T^g(fc-g)sinhfe 
PA,B,c,>{2kT,2qT,T) = ^^-^ (B.5) 

dz{k — qz)^ p /kz — qk — qz\ 
.1 cosh — cosh ' ' \k — qz^ 2q /' 

After the change of variables p = p(k) but before the change of variables x = x{z), the sub- 
threshold part, p^{2kT,2qT,T), has the same expression as /?>. However, in order to perform the 
change of variables x{z), the integral must first be split, J^-^ dx = ^^i^'^ dx+ J^^g dx. These integration 

intervals are mapped to the 2;-integrals and . Because ^ < on these two intervals, a minus 
sign relative to p> appears. One thus finds 

T*e(q-k)smhk / f-^ f'^\ , 
dz{k — qz)'^ [kz — q k — qz 



(kz-q k-qz\ 



cosh qz — cosh A; ' ' \k — qz'' 2q 

For the contribution 



2^5 ri ^ p cosh(p + ^/p"^ + + 2pqx)T — cosh(p — ^/p'^ + q'^ + 2pqx)T 

J -I Jo cosh{p + ^y p'^ + q'^ + 2pqx) — cosh.{p — ^y p"^ + q'^ + 2pqx) 



we follow exactly the same steps as for A, B, C and find 

, _ r4^(fc-g)sinhfe , 
f>v,> {2kT, 2qT, T) = ^-^^^^ (B.8) 



/. 



dz{k — qz)"^ {k"^ — q^ z^) /kz — q k — qz 



p^f l^z - q fc - qz \ 
\k — az' 2a / 



_i cosh A; — cosh \k — qz 2q 

For the sub-threshold part (which includes the minus sign in the numerator of Eq. B.7) we find the 
same expression as for p> before the change of variables x{z). Therefore, again, in the final expression 
/9< has a relative minus sign as well as the complementary integration range as compared to p>, 

. ,^ T^9(q - k) sinhk / f-^ f'^\ , 
pv.<(2kT.2,T,T) = + / ) (B.9) 

dz{k — qz)'^{k'^ — q^z"^) ^ /kz — q k — qz 
cosh qz — cosh k \k — qz^ 2q 

Finally, in view of the symmetry P{—y, —k) = P{y, k) of the relevant polynomials, /^^^ can be replaced 
by2/,^and/_-i + /rby2/-. 



C. The spectral function in terms of polylogarithms 



Polylogarithms are defined by the series Li„(2;) = Ylk>i W ^"^^ analytic continuation thereof. The 
integrals appearing in the spectral function are given as follows in terms of these special functions: 

q3T([z2];2rw,2rq,r) = -q^ ((log(l - e^"") - log(e-+'i - 1)) (C.l) 

-2q (Li2(e'i-") - Li2(e"+^)) + 2 (Ligle'i-*) - Li3(e-+^)) 

+iTTc^ -]- (1 + 6l(q - w)) (w^ - 27r^w + Szttw^) 
o 

q5X([/];2rw,2rq,r) = -q^ (log(l - e^-") - log(e<i+- - 1)) (C.2) 

-4q3 (Li2(e^-") - Li2(e"+'i)) + 12q2 (Li3(e^--) - Li3(e^+-)) 
-24q (Li4(e^-") - Li4(e^+™)) + 24 (Li5(e<i-^) - Li5(e<i+^)) 

- ^ (1 + 0(q - w)) (3w^ - 207rV^ - Stt^w + ISzttw'') 

The integral with in the numerator can be expressed in terms of elementary functions, see Eq. 2.36. 
Hence in all cases the expressions on the left and on the right of the threshold differ by a polynomial 
in 00. In order to simplify expressions Eq. (C.l) and (C.2), we have used the identities 

Li,(e--) - Li,(ew) = ^ _ - |lw + gw^ (C.3) 

Li3(e-™) - Li3(e-) = ^ - + y w^ . (C.4) 
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